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ABSTRACT 

We report the first observation of an anisotropy in the arrival direction of cosmic rays with energies 
in the muhi TeV region in the Southern sky using data from the IceCube detector. Between June 
2007 and March 2008, the partially-deployed IceCube detector was operated in a configuration with 
1320 digital optical sensors distributed over 22 strings at depths between 1450 and 2450 meters inside 
the Antarctic ice. IceCube is a neutrino detector, but the data are dominated by a large background 
of cosmic ray muons. Therefore, the background data are suitable for high-statistics studies of cosmic 
rays in the Southern sky. The data include 4.3 billion muons produced by downgoing cosmic ray 
interactions in the atmosphere; these events were reconstructed with a median angular resolution of 3 
degrees and a median energy of ~ 20 TeV. Their arrival direction distribution exhibits an anisotropy 
in right ascension with a first harmonic amplitude of (6.4 ± 0.2 stat. ± 0.8 syst.) x 10"**. 

Subject headings: cosmic rays — neutrinos 
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1. INTRODUCTION 

Long-term observations of cosmic ray muons by un- 
derground experiments have demonstrated the pres- 
ence of an anisotropy in the cosmic ray intensity up 
to a few hundred GeV (jNagashima et al.l I1998D . Re- 
cent underground and surface array me asureme nts of 
cosmic rays by the Tib et Array (|Am cnomori et aLl 
12006), Supcr-Kamiokandc ( Guihian et aH 2007) and Mi- 
lagro ( Abdo et al. 2009) indicate that the anisotropy per- 
sists into the TeV range. 

All of the TeV measurements were performed in the 
Northern hemisphere; so far, no such measurement has 
been performed covering the entire Southern hemisphere 
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at median energies in the multi TeV region. With the 
deployment of the IceCube Neutrino Observatory at the 
South Pole, we have for the first time measured the 
anisotropy at TeV energies in the Southern sky. Ice- 
Cube is primarily a neutrino detector but it is sensitive 
to the muons produced in downward-going cosmic ray air 
showers. The observatory provides high-statistics mea- 
surements of cosmic rays with median energy of 20 TeV. 
When completed in 2011, IceCube will comprise 5160 
optical modules buried 1450 and 2450 meters below the 
surface of the polar ice sheet. The modules are physically 
connected to the surface by electro nic umbilical lines, o r 
"strings," with 86 strings in total (jAbbasi et al.l 120091) . 
In this paper, we use cosmic ray data recorded by the 
detector in its 22-strings configuration (IC22) between 
June 2007 and March 2008 to produce the cosmic ray 
skymap of the Southern sky in the TeV range. 

2. ANALYSIS 

During the IC22 physics run, cosmic ray events were 
observed at an average trigger rate of about 550 Hz. The 
arrival direction is determined by a likelihood based re- 
construction which is se eded with a fast on line estimate 
of the arrival direction (jAhrens et al.|[20Q^ . The likeli- 
hood based reconstruction is applied if twelve or more 
optical sensors on at least three strings were triggered 
by the event. A total of 5.2 x 10^ events satisfied the 
above conditions at an average rate of ^ 240 Hz. Fur- 
ther selection criteria were applied to the data to ensure 
good quality and stable runs. The final data set contains 
4.3 X 10^ events with a total livetime of 226 days, a me- 
dian angular resolution of 3°, and a median energy per 
cosmic ray of 20 TeV. The energy scale was determined 
with a standard cosmic ray simulation program, COR- 
SIKA (|CORSIKAll2009l) . usingjhc SIBYLL hadronic in- 
teraction model (Version 2.1) (|Eii gcl 1999) and the Poly- 
Gonato model for th e composition a nd spectrum of the 
primary cosmic rays (|Horandelll2003[ ). 

To evaluate physical anisotropies in the cosmic ray 
data set, it is necessary to eliminate spurious effects 
which can mimic an anisotropy. These include local ef- 
fects such as diurnal and seasonal variations of atmo- 
spheric conditions, asymmetries in the detector geome- 
try, and nonuniform detector exposure to different re- 
gions of the sky. Fortunately, the location of IceCube 
at the South Pole is ideal to compensate for many ef- 
fects that can impact cosmic ray detectors in the mid- 
dle latitudes. At the South Pole, the Southern celestial 
sky is fully visible at any given time, providing com- 
plete and uniform coverage. While the seasonal vari- 
ation i n the cosmic ray event rate is on the order of 
±10% (jTilav et al.ll2010r ). these variations are sufficiently 
slow to have no effect on the anisotropy. Rapid atmo- 
spheric changes which can affect the rate are rare and 
can be identified from the data. 

The remaining effects which must be accounted for in 
this analysis are an asymmetry in the IceCube detector 
response, and a non-uniformity in the time coverage of 
the data. The asymmetric response is due to the geomet- 
rical configuration of IceCube during the IC22 physics 
run (as shown in Figure [ij; events arriving along the long 
axis of the detector were preferentially selected by the on- 
line filter and reconstruction due to the larger number of 
strings and modules triggered. In principle, the rotation 



of the Earth should average out the local asymmetry in 
the arrival directions each day, but gaps in the detector 
uptime and uneven run selection due to quality selection 
introduce non-uniformities into the time coverage of the 
data. These non-uniformities preclude the complete av- 
eraging, and translate into an artificial arrival direction 
asymmetry in equatorial coordinates. 
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Fig. 1. — This plot shows the IceCube detector geometry in the 
22 string configuration. The filled green circles are the positions of 
IceCube strings and the filled blue circles display the position of 
the IceTop tanks. 



To correct for this detector-related asymmetry, each 
event from a given local azimuth bin i was weighted 
with the ratio n/ui, where n is the average number of 
events over the full range of local azimuths, and Ui is 
the number of events in local azimuth bin i. Since the 
local azimuth distribution varies with zenith angle, the 
events were grouped into four zenith bands with approxi- 
mately equal numbers of events per band. The weighting 
is applied within each band to remove the detector asym- 
metry. 

3. RESULTS 

To investigate the arrival direction distribution of the 
cosmic rays, we studied the relative intensity of the cos- 
mic ray induced muon flux. The arrival direction distri- 
bution is dominated by the zenith angle dependence of 
the muon flux. The zenith angle dependence is a result of 
varying overburden for the muons through the Antarctic 
ice. Therefore, the flux was normalized within declina- 
tion belts of width 3°, which corresponds to the angular 
resolution of the data. This procedure provides the rel- 
ative intensity of the event rate in each declination belt 



independently. 

Figure [H shows the relative intensity of the event rate 
in equatorial coordinates. The color scale quantifies the 
number of reconstructed events with respect to the av- 
erage number of events in each declination belt. Fig- 
ure [3] shows the same data in Galactic coordinates. Note 
that since the declination belts in the equatorial map are 
treated independently, the map provides only informa- 
tion on the relative modulation of the arrival direction of 
cosmic rays along the right ascension. 

Figure [2] shows an anisotropy that appears to be a con- 
tinuation of a similar modulation of t he cosmic ray flux 
obseryed in the Northern hemisphere ( Amenomori et ahl 
[200l [Guiihan et al.l [20071: lAbdo et al.ll2009D . To auan- 
tify the scale of the anisotropy, we fitted the right ascen- 
sion dependence of the data to a first- and second-order 
harmonic function of the form 



n=2 



AiCOs{i{a - (j)i)) + B, 



(1) 



where {Ai,(j)i) are the amplitude and phase of the 
anisotropy, a is the right ascension, and i? is a con- 
stant. Figure |4] shows the anisotropy profile in right 
ascension obtained by accumulating the relative inten- 
sity distributions from the declination belts. The er- 
ror bars are derived by propagating the statistical er- 
rors from each declination belt, and the gray band in- 
dicates the estimated spread from the fit values of the 
stability tests. The solid line indicates the fit of equa- 
tion ^ to the data. The first- and second-harmonic 
fit parameters to the one-dimensional projection in Fig- 
ure H are Ai = (6.4 ± 0.2 stat. ± 0.8 syst.) x 10"^, 
6i = 66.4°±2.6° stat. ±3.8° syst., A2 = (2.1±0.3 stat.i 



0.5 syst.) X 10 and ( 
with xVdof= 22/19. 



-65.6°±4.0° stat. ±7.5° syst., 




Fig. 2. — The relative intensity of the cosmic ray fiux in equato- 
rial coordinates. 



To estimate the stability of the result and the cor- 
responding systematic uncertainties we performed two 
types of tests: we checked whether the observation is 
stable against the choice of the particular event sample 
selection, and whether the modulation has spurious in- 
fluences from other physical effects. The stability of the 
result was tested through a series of dedicated checks. 
The first stability test was done by dividing the data set 
in half by sub-run number, where each sub-run contains 
approximately 20 minutes of observations. To avoid any 
systematic biases the division was tried in several ways: 
first, by separating data in even- and odd- numbered sub- 
runs, and second, by random selection of half of the sub- 
runs. The corresponding relative intensity distributions 
in right ascension for both tests were determined, and 
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Fig. 3. — The relative intensity of the cosmic ray flux in Galactic 
coordinates. 
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Fig. 4. — The one-dimensional projection in right ascension 
of the two-dimensional cosmic ray map in equatorial coordinates. 
The data are shown with statistical uncertainties, and the black 
line corresponds to a fit to the data. The gray band indicates the 
estimated spread from the fit values of the stability tests (see text). 



it was found that the variation induced by the data set 
selection are within the statistical fluctuations. In addi- 
tion, to check for daily variational effects, the data were 
divided in two sets: the first containing sub-runs with 
event rates above the median value for the corresponding 
day, and the second containing sub-runs with event rate 
below the median value. The corresponding relative in- 
tensity modulations were fit and found to be smaller than 
the statistical fluctuations. This means that a variation 
in the absolute event rate does not affect the modulation 
on arrival direction in right ascension. 

Similarly, more stability tests were applied to check 
for effects due to seasonal variations and time gaps. To 
check for the seasonal effect data were divided into one 
set containing the winter months (June-October) and one 
set containing summer months (November-March). The 
relative intensity variations in right ascension were fit 
and found to be consistent with the statistical fluctu- 
ations. To verify that the non-uniform time coverage 
due to missing sub-runs and other gaps in the data is 
correctly handled by the azimuthal re-weighting proce- 
dure (see Section 2), the relative intensity distribution 
from the full data set was compared with the one deter- 
mined using only the days with minimal time gaps. The 
differences were found to be consistent with statistical 
fluctuations. 

In each of the above stability tests, an independent fit 
was made to the relative intensity distribution as a func- 



tion of right ascension using Equation [TJ The envelope 
from all the stability tests fit curves was constructed and 
it is shown as the gray band in Figure SI 

To verify whether the analysis procedure could in- 
duce a modulation in right ascension, the experimental 
event arrival directions were randomized to generate an 
isotropic distribution, and the same analysis was per- 
formed on this sample. The result was found to be con- 
sistent with isotropy. 

To check whether the observed anisotropy has some 
sidereal spurious effect derived from the interference be- 
tween possible yearly-modulated daily variations, the 
same analysis was performed using the anti-sidereal time 
frame (a non-physical time defined by switching the 
sign o f the transfo rmation from universal to sidereal 
time) (jFarlev et al.l 11954). The real feature in the side- 
real time is expected to be scrambled in the anti-sidereal 
time. Figure [5] shows the one dimensional projection in 
right ascension for the sidereal time in black and for the 
anti-sidereal time in red. The amplitude of the first har- 
monic fit to the one dimensional projection in the anti- 
sidcrcal time was found to be 0.8 x 10~^. This value is 
larger than the spread found in the first harmonic ampli- 
tude from the stability tests, therefore we use it as the 
systematic uncertainty in the first harmonic amplitude. 
The uncertainty in the first harmonic phase implied by 
the study in the anti-sidereal time frame is within the 
systematic error determined from the stability tests. The 
systematic uncertainties for the rest of the parameters of 
the fit (quoted in Sec. 3) are derived from the stability 
tests. 
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Fig. 5. — The one-dimensional projection in right ascension for 
the data in sidereal time (in black) and in anti-sidereal time (in 
red). The black line is the fit to the sidereal modulation using the 
first and second harmonic fit of Equation [T] The red line is the fit 
to the anti-sidereal modulation using the first harmonic only. 



4. DISCUSSION 

Using a high-statistics sample of downgoing cosmic 
rays with a median energy of 20 TeV and a median angu- 
lar resolution of 3° , we presented the first map of the rela- 
tive intensity of the flux of TeV cosmic rays in the South- 
ern sky. The arrival direction distribution of the cosmic 
rays is found to be anisotropic with a first harmonic am- 
plitude and phase of Ai — (6.4 ± 0.2 stat. ± 0.8 syst.) x 
lO"'^ and (/)i = 66.4° ± 2.6° stat. ± 3.8° syst.. The ob- 
servation appears to be a continuation of a previously- 



measured co smic ray anisotropy reported in the Northern 
hemisphere (lAmenomori et alJl2006tlGuillian et al.l[2007l : 
lAbdo et aLll2009l) ! 

The origin of the anisotropy remain s unclear. 
Compton-Getting (ICompton fc Gettind ll935') suggested 
that the relative motion of the solar system around the 
Galactic center in the cosmic ray plasma should give rise 
to an excess in the direction of motion of the solar system 
and a deficit in the opposite direction. In this model, an 
excess flux should appear with a maximum in right ascen- 
sion between 290° and 340° an d a minimum in right as - 
cension between 110° and 160° (jAmenomori et al.ll2006l ). 
As shown in Fig 4, the excess can not be described in 
terms of the direction of motion of the solar system. 
Therefore, we conclude that the Compton-Getting effect 
could be (at most) one of several contributions to the cos- 
mic ray anisotropy. This effect will be addressed in more 
detail in a future study including the energy dependence 
of the anisotropy. 

It is tempting to try to interpret the cosmic ray excess 
as an artifact of the heliospheric magnetic field. However, 
the maximum gyro-radius of a 10 TeV cosmic ray pro- 
ton in a 1 /zG magnetic field is about 0.01 pc, i.e. much 
larger than the size of the heliosphere. As a consequence, 
the observed anisotropy is more likely to be connected to 
features of the local interstellar magnetic field at dis- 
tances < 1 pc. We are also investigating the possibility 
that the cosmic ray excess is associated with structures 
in the Galactic magnetic field at larger distance scales, 
or with diffusive particle flows from a nearby Galactic 
source such as Vela. 

The still growing IceCube observatory will be com- 
pleted in 2011 with a total of 86 strings and a volume of 
1 km^ . The estimated rate of cosmic ray- induced muons 
will be greater than 30 billion events per year. Such high 



statistical power, together with an estimated energy res- 
olution of about 0.3 in log(E), will allow us to determine, 
in one year, the variation of cosmic ray anisotropy in sev- 
eral energy ranges up to a few hundred TeV. The energy 
dependence study will provide fundamental hint at the 
nature of the source or sources of the cosmic rays, as well 
as their propagation through the Galactic magnetic field. 



5. ACKNOWLEDGEMENTS 

We acknowledge the support from the following agen- 
cies: U.S. National Science Foundation-Office of Polar 
Program, U.S. National Science Foundation-Physics Di- 
vision, University of Wisconsin Alumni Research Foun- 
dation, U.S. Department of Energy, and National Energy 
Research Scientific Computing Center, the Louisiana 
Optical Network Initiative (LONI) grid computing re- 
sources; Swedish Research Council, Swedish Polar Re- 
search Secretariat, Swedish National Infrastructure for 
Computing (SNIC), and Knut and Alice Wallenberg 
Foundation, Sweden; German Ministry for Education 
and Research (BMBF), Deutsche Forschungsgemein- 
schaft (DEC), Research Department of Plasmas with 
Complex Interactions (Bochum), Germany; Fund for 
Scientific Research (FNRS-FWO), FWO Odysseus pro- 
gramme, Flanders Institute to encourage scientific and 
technological research in industry (IWT), Belgian Fed- 
eral Science Policy Office (Belspo); Marsden Fund, New 
Zealand; Japan Society for Promotion of Science (JSPS); 
the Swiss National Science Foundation (SNSE), Switzer- 
land; A. Kappes and A. Grofi acknowledge support by 
the EU Marie Curie OIF Program; J. P. Rodrigues ac- 
knowledge support by the Capes Foundation, Ministry 
of Education of Brazil. 



REFERENCES 



Abbasi, R., et al. 2009, Nucl. Instrum. Meth., A, 601, 294 
Abdo, A., et al. 2009, Astrophys. J., 698, 2121 
Ahrens, J., et al. 2004, Nucl. Instrum. Meth., A, 524, 169 
Amenomori, M., et al. 2006, Science, 314, 439 
Compton, A. H., & Getting, I. A. 1935, Physical Review, 47, 817 
CORSIKA. 2009, http://www-ik.fzk.de/corsika/ 
Engel, R. 1999, in International Cosmic Ray Conference, Vol. 1, 
415 



Farley, F., et al. 1954, in Physical Society. A., Vol. 67, 996 
Guillian, G., et al. 2007, Physical Review D, 75, 062003 
Horandel, J. R. 2003, Astroparticle Physics, 19, 193 
Nagashima, K., et al. 1998, Journal of Geophysical Research, 103, 

17429 
Tilav, S., et al. 2010, ArXiv;astro-ph/1001.0776 



